Brève introduction sur les points singuliers des équations différentielles et la théorie de 
Frôbenius pour les développements en série infinie autour de ces points singuliers. Éléments de 


Théorie Fuchsienne pour les équations différentielles du second degré 


Mon but ici est juste d'introduire quelques résultats de la théorie des équations différentielles 
présentant des points singuliers et de la façon de les résoudre en général. Cette théorie est aussi 
communément appelée Théorie Fuchsienne des équations différentielles du nom de son principal 
initiateur le Mathématicien Lazarus Immanuel Fuchs. Je me réfère essentiellement au cours de 
E.L.Ince dans son ouvrage « Ordinary Differential Equation a pages 356 et suivantes, page 396 et 
suivantes « Solution of linear differential equations in series », dédiés aux équations différentielles 
ordinaires d'ordre n. J'utilise également les notations introduites dans l'ouvrage « Field Theory 
Handbook » des auteurs P. Moon et E.Spencer dont un chapitre est dédié à la théorie de Frôbenius. 


Une équation différentielle linéaire du second degré peut le plus généralement s'écrire sous la 
forme : 


polz" ()+ pi (z)'(z2)+ 2,()y()=0 
Où les fonctions P? (z) 6) P-E) sont a priori quelconque. 


Les points « non-sinquliers » 


Les points singuliers se définissent en quelque sorte par opposition aux points ordinaires ou encore 
« non-singuliers ». Le point « non-singulier » est un point autour duquel la solution y(z) de 
l'équation différentielle peut être développée en série de Taylor, à partir de la seule connaissance 
de la fonction et de sa dérivée première en ce point : 


E Lutz, 

y(z)= x(z0)+ Le -zo )y'(z0)+ 

Autour de ce point, on peut alors construire deux solutions indépendantes, comme suit : 
ee (2)=1 EE 
»(z)= OCH (z SCHT Glt »() ES Y2(o)+ (z -=z )y' (2)+ Kä 


Autrement formulé c'est également le point autour duquel la solution du problème de Cauchy pour 
l'équation différentielle existe (problème aux valeurs initiales connues) : 


zm Yo 


lees gl 


Les points « singuliers » 


Par opposition les points « singuliers » sont ceux autour duquel on ne peut construire de série de 
Taylor pour la fonction y(z) (ou que l'on ne peut résoudre le problème de Cauchy). Pour caractériser 
ces derniers on peut supposer le plus généralement qu'il existe un développement autour du point 
singulier d'une forme particulière. On parle également d'une fonction solution de l'ED d'ordre fini 
autour du point singulier lorsque qu'il existe une valeur p finie telle que la limite suivante (pris dans 
le plan complexe) est nulle ` Lim (z-2z,) y(z)=0. Pour le point à l'infini, cette dernière est définie 
. . Lim z°y(z)}=0 — , ss 7 x EN 

ainsi: =>> . Si les deux solutions indépendantes sont réputées d'ordre fini autour du 
point singulier, alors le point singulier est dit régulier et irrégulier si ce n'est pas le cas. 


Voyons comme cela se concrétise avec l'équation différentielle définie plus haut : les seuls points 
singuliers sont ceux pour lesquels p,(z) s'annule dans le plan complexe, à l'exception du point à 
l'infini. Pour simplifier le propos, il est d'usage de placer un de ces points singuliers en z=0, sans 
pour cela réduire la généralité du propos. Le développement est donc supposé de la forme : 


SE 2 
tel ser Sek Dans ce cas il suffit de prendre p=1-r pour que la limite précédente 


s'annule et que la fonction soit réputée d'ordre finie. 


Les fonctions Po (z) AG) ?:G) sont également supposées être des fonctions analytiques dans tout 
le plan complexe, y compris aux abords du point singuliers z. Elles admettent donc chacune des 
développements de Taylor. En injectant les développements propre des fonctions 20770770) 
per) i 2 e D D e D 

et le développement * ()=7 (C +Cz+0C2° +...) dans l'équation différentielle, et en annulant 
chacun des termes de puissance de z, on obtient un système infini d'équation linéaire sur les 
coefficients C; qui s'annule. La première équation linéaire de la puissance minimale en z est appelée 
communément l'équation indicielle. Dans le cas d'une équation différentielle du second degré c'est 
un polynôme de degré maximal 2, que je note I(r). 


Plusieurs cas permettent une première classification du point singulier z=0 : 


Cas n°1 : I(r) est indépendant de l'exposant r, il n'y a pas de solution et le point est alors qualifié 
d'irrégulier 


Cas n°2: I(r) dépend de l'exposant r, si c'est un polynôme de degré strictement inférieur 2, là 
encore le point z=0 est un point singulier irrégulier 


Cas n°3 : I(r) dépend de l'exposant r, c'est un polynôme de degré = 2, le point z=0 est un point 
singulier régulier. 


Pour le cas n°3, modifions l'équation différentielle de la façon suivante : 
y"(z)+ P(z)y'(z)+ Q(z)x(z)=0 P()- PE) och 20). Pour qu'un point singulier z=0 de 
Po(z) Po(2) 
l'équation différentielle se trouve dans le cas n°3, il faut et il suffit que les deux fonctions P(z) et 

Q(z) satisfassent aux conditions suivantes (Théoème dû à Fuchs) : 


P(z)= 22 +0() 


EE ET 


2 e 
Alors dans le cas n°3 l'équation indicielle est une équation du second degré 7 +r(pe—1)+4 = 0, 


Elle s'obtient facilement en injectant les développements approximées de P(z) et Q(z) et en 
annulant le coefficient d'ordre minimum en z, si l'on suppose un développement commençant par 
le terme z'. On voit donc que l'équation indicielle s'obtient à partir du calcul des limites suivantes : 


; ; PE) 

=L P\z)=L A 
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=li 2 = Li 
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Définie comme cela, on peut encore reformuler ces limites comme un calcul de résidus de fonction 
de valeur complexe pour un point singulier régulier quelconque 2, : 


Po = Residu —— | 


qo = Residu … (e-a) 2) 


L'équation indicielle admet deux racines ro et rı que l'on appelle communément exposants 
caractéristiques de Frôbenius. Il y a alors trois cas possibles dont nous allons donner les résultats 
principaux. Par la suite les développements seront explicitement construits autour d'un point 
singulier z0, à l'image du point singulier z=0. 


Sous cas n°1 : les deux racines ro et rı sont distinctes et l'écart entre les deux racines n'est pas un 
entiers, alors on peut construire deux solutions indépendantes autour de z=0, sous la forme : 


l=+%0 


1=+0 
rer Cz wkl? > Dz 
1=0 1=0 


Sous cas n°2 : les deux racines ro et rı sont distinctes et l'écart entre les deux racines est un entier, 
alors on peut construire deux solutions indépendantes autour de z=0, sous la forme : 


1=+0 Jeton 
Yy(z)=z* KS »(2)= »,(z)Log(z)+z° 3 D,z' 
1=0 1=0 


On observe la présence d'un terme dit logarithmique. 
Sous cas n°3 : les deux racines ro et rı sont identiques, alors on peut construire deux solutions 


indépendantes autour de z=0, également avec un terme logarithmique : 


1=+00 1=+o 
Her Cz ()= yo(z)Log(z)+z* 3 Dei 
1=0 


1=0 


Le points « singulier » à l'infini 
Il est facile de transformer le point à l'infini en un point singulier z=0 alors la transformation z->1/z. 


L'équation différentielle de départ devient par cette transformation l'équation différentielle 
suivante : 


y"(z)+ PG)r'G)+0(G)y()= 0 
c=1  y()=(ç) 


Z 


ant. Allan. Eu. 


Les conditions pour lesquelles le point à l'infini est un point non-singulier sont alors les suivantes : 


Les conditions pour lesquelles le point à l'infini est un point singulier « régulier » sont alors les 


suivantes : 
EEN 
S ç (a (a G 


Dans ce dernier cas l'équation indicielle a la forme suivante, à partir d'un développement des 


fonctions P(1/0) et OO Al : 


> r?+r(2- p, -1)+4 =r?+r(1- Date = 0 


Moyennant cette transformation le point singulier à l'infini est alors traité comme tout autre point 
singulier. Les valeurs pO et q0 se déduisent aussi d'un calcul de résidu : 


Po = Residu _, (28) 


Po) 


qo = Residu _., , Gel 


Po(z) 


Exemples d'équations indicielles dans le cas de points ordinaires, singuliers réguliers et singuliers 
irréguliers 


Supposons que dans l'équation différentielle le point z=z, soit un point ordinaire, alors on peut 
écrire : 
»"(2)+P6h'()+06)y()=0 
= Social" ` th leisst" yo- Nr sa? 
P(z)= Palrls ae lëlalt—  Q()=0()+(2-2)2" (20) 


Puissance minimale en z terme y"(z) — (r + I(r KH 1) 1=0 


Dans ces conditions l'équation indicielle est invariablement 7 (@-1)=0 pour un point ordinaire ou 
non-singulier avec deux solutions r=0 et r=1. 


P(z) z= 


Po_+o(1) 
On a vu que pour un point singulier telle que : ge , l'équation indicielle était 


1 


— Zo 


o(z)= G Se SP 


2 pes 2 . 
de la forme :" + (Po —1)+4 =0 Le point z=z, est un point singulier régulier. 


P(z)= Po 0] l 
z=z) Z—2Z9 


Supposons maintenant que d alors exprimons juste le polynôme indiciel 


QE ee LE | 


I(r) du degré minimal en puissance de z-z,: 
eOVOE GC St 4) qo 
P(z)y'(z) =r C, (z =i y> Po — Puissance minimale en z terme P(z)y'(z) — r=0 


y"(z) —> r(r — 1Xz — Ze des Ci 


I(r) est un polynôme de degré strictement inférieur à 2. Donc z=Z9 est un point singulier irrégulier. 


P(z)= Po o cl l =) avec n>l 
Pour une singularité de la forme : (2-20) C-z) , alors polynôme 


o-de] E 


indiciel I(r) est de la forme : 


Qlz)y(z)= C (z- z0) "o Si n+l<m > I(r)=% 
P(z)y'(z) =r C; (z SZ; e Po — Puissance minimale enz terme Si n+l=m — I(r) =F Po + qo 
y"(z)—>r(r-1Xz-z)°C, Si n+l>m >I(r)=r 


Là encore le z=z, est un point singulier irrégulier. 


Notation de Riemann pour les points singuliers réguliers d'une équation différentielle du second 


degré 


Lors de l'étude de l'équation différentielle suivante (aussi appelé équation de Riemann ou de 
Papperitz) : 


+ 
z-a z—b Z—C 


UE ET 1-B-f 1-7-7 a œ'(a—-b\a-c '(b—-ab-c '(c—afc—-b 
vil LIB=B" 1-7 "kd (a-b) log X-c), r'(c-aX d 
z—a z—b ET (z-az-b\z-c) 
Une notation commode a été introduite pour désigner les solutions formelles de cette équation. Ici 
les points z=a, z=b et z=c sont des points singuliers réguliers d'exposants respectifs pour z=a, a,a' 
pour z=b, 6,6' et pour z=c, y, y'. La solution est alors notée symboliquement sous la forme : 


dénommé symbole P-Riemann ou symbole de Papperitz (voir la combinaison des deux noms). 


Les trois premières colonnes donnent respectivement le point singulier régulier et les deux 
exposants solutions de l'équation indicielle lors du développement autour de chaque point 
singulier. Par extension ce symbole est particulièrement adapté pour représenter formellement la 
solution d'une équation différentielle à nombre quelconque de points singuliers réguliers, sans se 
soucier de sa résolution concrète. On peut même lui appliquer le point singulier à l'infini (s'il est 
régulier) et une transformation de l'argument z. 


Pour l'exemple prenons l'équation différentielle de Bessel : 


Prenons un autre exemple avec l'équation de Heun (j'applique également la relation de Fuchs pour 
déterminer l'équation indicielle au point à l'infini): 


ebe ae ich ea pour 


z z-l z-a z-l{z-a 


Point singulier z=0 P=% qo 0>(r)=r+r(y-1)>r=0 r=l-y 


Point singulier z=1 p= qo 0= 1(r)=r° +r(8-1)= r =0 r=1-6 


Point singulier z=a Du zs E Qo 0= 1(r)=r?+r(e-1)= r=0 r=l-E 


Po = Lim EI à Le }rr+o+si-n--a-8 


zm z z-l z-a 


HD SE Le 
To Pe f z(z-1X{z-a) SS 


Point singulier z= = r?+r(1-p,)+q,=(r-ar-B)=0=r-a r=ß 


0 1 a œ 
On a donc le p-Symbole de Riemann : 4-1. P| 0 0 O a z 


l-y 1-6 l-e $ 


Équation « fuchsienne » et relation de Fuchs 


Une équation différentielle ordinaire d'ordre quelconque est dies fuchsienne si tous ces points 
singuliers, y compris le point à l'infini sont réguliers. 


Il existe une condition pour que tous les points soit réguliers. Je la donne pour une équation 
différentielle d'ordre 2 (une démonstration partielle simple suit). Soit a; l'ensemble des v points 
singuliers de l'équation à distance fini (en dehors du point à l'infini), soit a;,a: les exposants aux 
points a; et 6,6- au point à l'infini, la relation suivante a lieu : 


D CA + a)+ B, + B, =v-1 Equation différentielle de degré n SS Q; + 5 B; = gab 


i=l il j=l j=1 


Prenons par exemple l'équation de Heun, alors la relation s'écrit : 


r=3-—1-y+1-6+1-E5+a+B=-2Sl+a+pB=7+0+e 


Cette relation se démontre « facilement a si l'on admet le théorème du même mathématicien 
Fuchs que l'on a déjà énoncé et qui nous dit qu'une condition nécessaire et suffisante pour qu'un 
point singulier soit régulier pour l'équation ;"(:)+ P(z)y'(z)+ QO(z)y(z)=0 est que les deux fonctions 


P(z) et Q(z) soient de la forme P(z)= PIE? Q(z)= CF autour du point singulier z, avec des 
Z= Zo Z— Zo 
fonctions po(z) et qo(z) analytiques (admettant un développement de Taylor) autour du point 


singulier. 


Si maintenant on transpose ce résultat dans une équation Fuchsienne, cela implique que dans la 
partie finie du plan complexe avec pour seuls points singuliers a;, alors les fonctions P(z) et Q(z) 
sont nécessairement de la forme : 


BES Pu +q(z) och Poi a) 


+ 
SE (z-a) G-a)x.x(-a) 


= 
Il 


Avec qm(z) des polynômes en z de degré à définir. 


Pour cela étudions maintenant le point singulier z=°. Pour que ce point soit également un point 
régulier, posons le développement de la solution à l'infini sous la forme: 
cl zle taz" res" +...) Et supposons encore que le développement à l'infini des fonctions P(z) 
et Q(z) soit de la forme : oi (bo+biz +.)  Qle)= 27hb +b,” +.) Calculons alors le 
polynôme indicielle de telle manière qu'il soit exactement du second degré (condition de régularité 
du point). Dans l'équation différentielle ,"(;)+ P(z)y'(z)+O(z)y(z)=0; les termes de puissance 
minimale en 1/2, soit les termes dominants en puissance de z sont: 

y"(z) >r (r = 1) 7e 
P(z)y'(z)> p z 
Ochizl-zsz 


Pour que l'équation indicielle soit du second degré il faut et il suffit que le puissance maximale 
vienne du terme de la dérivée seconde. On a donc les conditions : 


En >r—-l+0, <r Ge m = {1,2} 


Q(G)r(z)- r+0, <r—2 


Cela signifie que finalement les polynômes qm(z) sont d'un degré inférieur ou égal à m(v-1)-v, soit 
q ()=o("t). Cela implique immédiatement que pour m=1 si qm(z) doit rester un polynôme 


alors a(z)= SIS et pour m=2, a ()=0("?) Si par exemple il y a deux points singuliers à 


distance fini, alors ce polynôme est une constante y 3, (:)=o(1)=a: Avec trois points 


singuliers alors, =3 = q,(z)=0(z)=a+b z: 


On en déduit la forme la plus générale d'une équation fuchsienne avec v points singuliers réguliers 
à distance finie : 


oga oir Sn 


l=1 ig SE 

Il reste à tirer la conséquence de la régularité du point à l'infini. Maintenant que tous les points 
singuliers à distance finis sont réguliers, les exposants caractéristiques de tous ces points avec la 
forme obtenue de l'équation fuchsienne sont solutions d'une équation du second degré de la forme 


Šp ER . 
cker put pau =r +r (pu 1)+ pa, et peuvent être pré-assignés aux valeurs : 1% (eau 


Mais comme ce sont les racines de l'équation du second degré, on peut écrire : 


Gs A > 
r? +r (pu Ap, =r- a)r aiz), On a donc les relations évidentes : 


l ZP = da + 


Pai 7 An 


L'équation différentielle au point à l'infini a la forme asymptotique suivante : 


i=v 


Pii 


BEE eich AB tan Jatci-g 


Partant de là l'équation indicielle au point singulier z=% a la forme suivante : 


deer Ba, Bn, H T H B,Xr H B,)=0 


Pi + B; =$ p; =] 
i=l 


On en déduit les relations suivantes : portant sur les exposants 6, et 62 à 


BB; = D Pa: TO 


i=l 


l- p; = +; 


Pai 7 hd 


l'infini. Mise en correspondance avec les relations précédentes d 2 , on déduit : 


< 
i 
< 


i=v i= 


(1 ol Lo, to lech E 1-8 b= X (ar +an)+ b +b, =v-1 
ai 


i= 


i= 


Ce qui est la relation de Fuchs recherchée. 


Au passage la seconde relation nous dit que pp, = Yaaa tha = Qu = BB - D Ana 
i=1l i=1 
Il est à noter que les coefficients 4, _,...,4,,,q2, PE Sont pas déterminés. Avec ces relations la forme 


de l'équation différentielle devient : 


EN a a 

S1l-a,-a, = er S > Se eg PULL 

" i > > S =0 
{ Ek HT Go) aad ea) ESCHER EEN us 


DS z-a a) 


CH La Lis v- v- 
z)+ "Dn KH Lia H (8,8, Ce +Q2-37 TA RÉCK 40) Far el 
a ZO 274; (z ol 
3 
Or (z Ce H TI a,)= Ta ol H (z as 32" ++ Liz + X20) 
IERCH IRCH 


Or cette forme peut encore se simplifier sachant que : 


i=v i=v 


(z Ada Il a,)= II te a,)+(z a (Zass? ++ X112 + X201) 


=li i=l,i#l 


où les coefficients y Lt, n'ont pas besoin d'être déterminé explicitement, mais dont on 


2-30 


sait que l'existence est assurée, car l'expression du membre gauche de l'égalité est égale à 


Ta, al lorsque z=a,. Donc la différence est nécessairement un polynôme qui s'annule pour z=a; 


ilit 


et qui est de degré égal à v-1. Il a donc nécessaire la forme : (z-a Zaraz" ++ aZ + Saach Une 


fois que l'on injecte cette forme dans l'équation différentielle fuchsienne, il vient : 


` oo IERCH 4 
rodga) yo 2 "Së xe) = 


l=v 


(z-a) 


=1 


l=v 


v-3 v-2 v-3 
+5 a ar 312 Geet GIE toi) (8,8, Z +537 ++ Qt) 


I=1 


l=v l=v l=v 
Posons Oe an + à He) = 21 = it ET EE 20 = 20 + D E 
E] pe CH 
i=v 
Sl-a,-a [Te -«) »G) 
il 12 i=li#l v-2 > v-3 ~ ~ Q 
> y"(z)+ > JS Outil ct (8,8, Z Fa +.. +424 Go) Se 0 
l=1 z-a Z= Gsa) 
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l=1 


Donc pour une équation fuchsienne où Pa; =0 une forme assez générale possible est la suivante : 


l=v v—2 
y"(z) | D Pi Lu , EE EES BH A 
` 


SES DEES 


On déduit facilement que : a, =0 o. 1-p,, €t également que Qu = BB Mais on va voir que 


la forme la plus « condensée » possible est obtenue à l'aide d'une suite de changement de 
fonctions. 


Si nous réalisons le changement de fonction suivant où m est à valeur quelconque réelle, il vient : 


J 
y ka: (z)+ 2m viel m(m 1) wc) -| 
zma, z-a, ) 
wt TTE 
wj g — i| An ell 17 di EH 7 7 8 
DND H b = -taye H >, S S H (88 SS TO +. +d Zt Lo) ren SS vu 
l=1 Tu l=1 "Hi za 
Il 
> w"(z): l y Lage, cerner) w'(z)+ 
IL j z-a CH, 
E = = Ch = l-a,-@; 
iy | And [Ila -a,) Q aQ j TIl, oh m(m D [IE a,) ; moll a,) X E 
Ziel A IERCH IERCH jst l=1 l í 
+| iej Z-a, | z-a; g el =0 


v-2 >~ v-3 pd pod 
BIZ? Z +4,37 en) 


> w"(z): l 5 Lage, "hel w'(z)+ 


ej ZO z-a; 


i=Lisl : =1,i# j i=l,iż j i=l Lisi 


iv i=v i=v i=v fey 1 Z SS i=v 
y | nänn [Ila -a,) Ga: TIl, a.) H m(m OU ol +m (fie ol SH On — y h o, Al il 
i=l i i= 


+ List Z-a z—a 
Mel. wë v-3 Gi Gi 
BI? Z "Fa enn) 


L'expression polynomiale suivante : 


Rats i (a, A WË § aan sha, Afen 


CH i=l,i# j 


l=v l=v 
Terme 2" uch (m-1)+m l-a, -a „+m 5 (l-a, a) 27 li -m+mÿ (l-a, all SE 


Lis? = 
P la,)= E +m(m-1)+m l-a, -a;) WÉI -a,)= lan "wl, alle a) 
fa Lie? 
> P (z)-P la )= 6-2 nn) (a, nl ge r = e-a fr -m+mv SmS Gs va) SC? 8 
= CG 
i=v l=v 


En utilisant la relation de Fuchs DC +0,)+8,+8B,=v-1= m°-m+v -mY (a; +a,)=m +m(B, + B,) 
i=l Wel 


SS P()-P{a,)- (e-a or +m(B, ls"? +.) 


Il vient une équation différentielle de la forme : 


> w"(z): | > bn Fe e, m) len ") w'(z)- 


Lisi Z7 Z—d; 


iy | Ann [Ila -a,) (a, -ma all Ia, -a,) 8 d e ( ) 
A? EES ' E l ID më, +m) z”? +4,32" Tt 3) — EH 


Bes z—a z-a, 
Ie j I j (z-a,) 


Donc la transformation y(z)=(z-a,}" w2) conduit à ce que l'équation est toujours une équation 
fuchsienne à v points réguliers à distance finie, pour lequel l'exposant au point z=a;, les exposants 


sont transformées de la manière suivante Ju" Ta ere, 
Q jz = Aam B:=B,+m 
En prenant tout simplement m=a;:, alors n T et l'équation différentielle devient : 


ORe 


DEE d ia aE- 
l 


Lisi Zoa Z— Ii z-a (z E ) 
1 


f= 


Donc l'équation différentielle c'est simplifiée en complexité notamment sur le terme : 


(age z-a 


En itérant ce processus pour tous les indices j entre 1 et v, soit en se fixant une transformation de 
IEN 


la forme : y(z)=(z-a,}" (z-a, (z-a, )" w(z)= well at: l'équation différentielle de la fonction w(z) 


l=1 


devient : 


vo l-a, +] ZG + Eau | KH SE E ON E wz) =0 


Se (z-a) 


Exposant Gu =0 Ost — An 


On peut donc considérer que la forme la plus standard de l'équation fuchsienne est la suivante avec 
le symbole de P-Riemann suivant : 


POSE) vas aa ttar 


Fa Ga 
I=1 
a a, a, œ 
P 0 o 0 Bz 


Une question pourquoi la forme la plus générale d'une équation fuchsienne à trois points 
singuliers est celle de Heun 


Appliquons ce résultat à l'équation à 3 points singuliers réguliers à distance finie : 


vof E Suite 


P 0 0 0 Bz 


l-a, l-a, 1-a, f, 


Il se trouve que cette équation différentielle reste de même nature par une transformation 
homographique simple z->z-a, (translation) : 


vil 3 E SH 


t=z-a = w"{t)+ Si SC } Fe w'(z)+ BiBat +a BB, —4 w(z)= 
ER k z) E ET 


Cette translation remplace le triplet de point (a ,a,,a,)>(0,a,-a,,a,-a,) et l'équation reste une 


équation différentielle de même nature après re-écriture des différents paramètres : 


SCHER Er Joch SE HE 


zZz z-a z-b DEET EE 


La deuxième transformation homographique est également très simple (homothétie) et associe les 
points Oas) [012] : 
a 


P EAS EEN E Jo BBSt-4 vn 
a orl ) 


at at-a at-b at-aYat-b 


Et l'on obtient également une équation différentielle de Heun « standardisée » sous la forme : 


EECH E te 


z z-l z-a z(z2-1{z-a) 


A propos d'une solution polynomiale de degré 1 de l'équation différentielle de Heun 


Soit l'équation différentielle fuchsienne comportant 4 points réguliers : 


vil e Een AS) 4-0 


Se Sch ig z(z-1X{z-a 
Contrainte de Fuchs a+fB+1l=7+ô0+e 


Je reviens sur la construction des solutions polynomiales de cette équation, en appliquant les 
résultat établit en 1931 par Moris Marden dans son article « On Stieltjes Polynomials ». Il s'agit de 
caractériser ou de construire les solutions polynomiales des équations fuchsiennes qui sont 
communément appelé « Polynômes de Stieltjes » ou « Polynômes de Stieltjes-Heine », ces derniers 
sont étroitement associés à l'existence et la forme de polynôme de Van-Vleck, présent en 
numérateur du terme en y(2) de l'équation différentielle fuchsienne (voir ci-après). Le paramètre n 
est le degré du polynôme de Stieltjes et p est le nombre de points singuliers réguliers à distance 
finie de l'équation fuchsienne. 


Et je me restreint tout d'abord à la solution polynomiale de degré 1, soit la plus simple qui se puisse 
trouver. Il se trouve que l'article laisse au lecteur le soin de compléter pour le cas n=1. Plaçons-nous 
également sur le cas p=3 soit l'équation de Heun. La forme du polynôme de Stieltjes dans ce cas 
peut s'écrire : S(z)=z-z En injectant cette forme triviale dans l'équation de Heun, il vient : 


y, 8 „£ , (@Bz-8) 


z z—-l z-a z(z-1\z-a) 


E OE a E EE L F di: e ) 


z zl z-a 


EE 


L'un des résultats de l'article est que par ailleurs la valeur de la racine du polynôme de Stieltjes est 
solution de l'équation algébrique de degré 2, qui est également une équation dites fractionnaire : 


y ô E 


H 
Z Z—-l z-a 


Ainsi l'expression z(z—1)\(z di GE ) s'annule en zën, Comme c'est une expression 
z z-]l z-a 


polynomiale de degré 2, on peut la factoriser sous la forme : 


z(z-1%z af 2 ER DE 


Il vient donc pour la détermination du polynôme de Van-Vleck : V(z)=a B 2-0 la valeur suivante : 


hr 0-2): e d SE E )- GE E EE 
z z—]l z-a 
= a Bz-06,=-(y+8+eXz-8B) 
a B=-(7+8+e)=(1+aX1+8)=0=ax=-1 où B=-1 
u ege ME 
y+ô+E l+a+ß a 
= V(z)=-B z-0, ou V(z)=-a z-0, 


Racine de V(t )=0 | poch ou E 
B a 


Le second résultat de l'article porte sur la valeur de la racine du polynôme de Van-Vleck qui pour 
n=1 doit suivre l'équation algébrique : 


L'expression de la racine z, est la suivante : 


fa Re re E 


H H 
DENGE 


ez (y+8+6) z(y(1+a)+8 a+e)+ay =0 


Lt e Iä-Zrebrzäbdn ,,-_(8-8+a(r+6)-4a 
2 B-5+a(y+6) -4a78 _ |P" 2 a 2 
A=(B-6 ô)a) -4 Ô +e)= 
ER E f +aļly +8) ger" An MÉ S+a(y+6))}+VA 3 __(a-6+a(y+6)-VA 
1 2 Z 2 
1_Bz-(8-8+(vr+6)a) 1 _-(8-6+a(y+6)+a_ po 
Pz -z(y+e+(y+6)a)+ay=0 | 2 ay 2. Ze ay "ay 
> > 
az -z(y+e+(y+6)a)+ay =0 1 _az,-(a-8+(7+6)a) 1 _-(a-8+a(y 6))+ VA ` Tue d 
z ay z 2ay ay on 


J'aurais pu encore tout aussi simplement observer pour les deux racines du polynôme de Stieltjes 
que: 


ylz ls, -a)+6 zl ele za -1)=0 (7+8 +e) zk a )=0>(7+8+e)z' z =ay 


; ay 

ZA ZA — D 
> 

y ay 

E 

a 


Soit avec la normalisation choisie pour les polynômes de Heun : 


(B-5+a(y+6))}+VA EN 


1+z— S =]+z—— ou l+z— 
Si (z)=z-z = siz); = ay mn Sé SS 
Z, a PR Ce gef, aa 
2ay ay ay 
a. Lë-ärebrzéëlizdä , ,-__(8-5+a(r+6)-Va 
i 2 ; 2 
avec 
+ (x-5+a(y+6))+VA = (xa—8+a(y+6)) VA 
d S et 0 = 2 


C'est évidemment l'expression des polynômes de Heun de degré 1. Ainsi pour ce cas particulier où 
n=1 et p=3 l'article de Moris Morden de 1931 donne bien comme polynôme de Stieltjes un 


polynôme de Heun. 


A propos des solutions polynomiales de l'équation différentielle fuchsienne générale 
Polynômes de Heine-Steiltjes et polynômes Van-Vleck 


Soit l'équation différentielle fuchsienne comportant p points réguliers : 


D 4 O SE EE 


LZ Hi 
1=1 


a, DH les deux exposants du développement de Frôbenius à z = 


l=p 
Contrainte de Fuchs a+fB+1= > Vi 


l=1 


V(z)= a B z”? tit? +... + Vo 


On va supposer que les coefficients y, sont tous réels et positifs. Dans l'article de 1931 de 
M.Marden « On Stieltjes Polynomials » la supposition sur les valeurs des coefficients y, s'étend à 
des valeurs complexes : _, < argy, <v < - , Soit le demi-plan complexe à droite, ce qui exclut les 


valeurs négatives. Un théorème dû à Heine (Handbuch der Kugelfunctionen Berlin édition de 1878, 

section 135, « Ueber die Existenz und Anzahl der Lamé’schen Functionen hoherer Ordnung », page 

473) nous dit qu'il existe au plus çr- = k SS _nx(n+1)x..x(r+p 2) polynôme V(2) de degré 
E p-2 1x2x...x(p-2) 

p-2 pour lesquels il existe des solutions polynomiales S(z) de degré n. S(z) est appelé polynôme de 

Steiltjes-Hein ou Steiltjes tout simplement et V(z) les polynômes de Van-Vleck. 


Un article de T.J.Stieltjes de 1895 « Sur certains polynômes qui vérifient une équation différentielle 
linéaire » dans « Acta Mathematica » page 41 et suivantes, l'auteur démontre qu'avec les 
conditions que nous avons choisi y >0 1-1...,p alors les déterminations des E polynôme 


V{2) de degré p-2 sont toutes réelles ainsi que les solutions polynomiales de Stieltjes. 


Concrètement lorsque n=1, il existe au plus p-1 polynôme de Van-Vleck de degré p-2. Et encore plus 
concrètement pour p=3, il en existe 2. C'est d'ailleurs ce que l'on a construit dans l'exemple simple 
précédemment exposé avec les deux racines possibles d'une équation du second degré. 


Le théorème de Heine établit l'existence de ces divers polynômes. Pour la construction de ces 
derniers, l'article de M.Marden en établit les prémices. Dans cet article le premier résultat porte sur 
les zéros du polynôme de Stieltjes de degré n. Soit les zéros zę coïncident avec l'un des points 
réguliers a,, soit ils satisfont à l'équation « fractionnaire » suivante : 


Pour n=1 k={l,. n> 0} 57 20 


11 Z Ou 


Le cas n=1 n'est pas traité dans l'article de M.Marden du fait de sa relative simplicité. De l'équation 
fractionnaire il en résulte par une mise en facteur commun un polynôme de degré p-1 au 
numérateur dont l'annulation détermine p-1 valeurs possibles de la racine z1. 


Pour déterminer les valeurs des polynômes de Vlan-Vleck, la construction est la suivante. Une fois 
déterminé les racines zx du polynôme de Stieltjes, alors on peut calculer les racines 2! de la dérivée 


première du polynôme de Stieltjes. Notons t, les zéros du polynôme de Van-Vleck de degré p-2 alors 
ces zéros soit coincident avec l'un des points réguliers o, soit ils satisfont à l'équation fractionnaire 
suivante : 


= Ki e 
Ki Pour n>1 k Dn 2} > 0 


Ee 
V(z)=a 8 IIG-:)- Se 4 "fil 
D Pour n=1 k= d'B 2} N His c 0 


Là encore le cas n= 1 étant considéré comme relativement simple par M.Marden il n'est pas traité 
et il résulte que les zéros t, suivent exactement la même équation fractionnaire que le zéro zı. De 
plus toutes les racines t, à choisir étant distinct, le problème revient à choisir p-2 racines distinctes 
de zı parmi ces p-1 valeurs, soit le coefficient binomial qui est justement égal à p-1. Il s'en suit qu'il 
y a bien p-1 formes possibles du polynôme de Van-Vleck V,(z). 


Une considération générale sur l'équation fuchsienne conduit à une condition nécessaire pour 


l'existence de solutions polynomiales 


Avant de passer au cas n=1 et n=2, une considération générale sur l'équation différentielle permet 
de déduire les seules valeurs admissibles de o. pour qu'une solution polynomiale S,(z) existe avec 
un polynôme de Van-Vleck de degré p-2 : 


S. (z)= (z z) l=p V (z) l=p l=p j=p 
E Projo sos ofie- nfe- s7030- 
V,(2)=a B [[G-1) pa (z-a) R FC 


(=D) la ener". sF SI lala Mer. 


1= 
Terme z”? = aln Däi +a B 
l=1 


l=p 


>n(n Ins Ee B=n(n 1)+n(1+a+8)+a B={(n+a]n+p)=0 n=-a ou n=-f 


Ern 
Relation de Fuchs KS =1+a+p 


I=1 


Polynômes de Heine-Steilties de degré 1 et Polynômes de Van-Vleck associés pour une «équation 


fuchsienne générale » 


Précisons donc la construction du polynôme de Heine-Stieltjes de degré 1. Nous avons illustré ce 
cas simple pour l'équation de Heun dans le paragraphe précédent. À titre d'exemple traitons ce 
même cas pour l'équation fuchsienne générale à p points singuliers réguliers : 


l=p 2 Lnl j=p Izp] j=p l=p l=p-1 

y(z)=S.(z)=z Zi 2e Ss =0 D [It -a,) =0—p-1 valeurs de AA la Jh DEEN 
= 474 = j=l = j=l = 1 

pél Gel 


l=p-2 


V(z)=a 6z”? +4 vpz” HtA JIC t)>4=aß 


l=1 


l=p-1 


p-1 formes deV(z) Hess) e)=a B IG k=1,..,p-1 


a p (z-z z-z) 
SE: $ l= l=p l=1 
TOR ge k=1, y% l A El a { Z Yı | (2): 0 ( ) 0 1 fi ï — Jet =0 
GH 1274 (z-a,) 1274 (z-a,) 
Ei LA 
l=p j=p 
TI -a) TI 
A ja ZS ) Lal rn I=p-1 
© = +a p— E =0 » 2 (z-a,) +a ß JIk z,)=0 
(z-a,) (z -a) Se Zo = 
l=1 I=1 
l=p l=p-1 l=p-1 l=p-1 


a, P les deux exposantsàz=% ot zs (a+B+1) (z-zu)+a p [IE 2,,)=0= (a -1X8 +1) 


(z z,)=>a= l1 ou B=-1 
l=1 l=1 l=1 I=1 


Donc il suffit que à ou 6 soit égal à -1 pour que le couple suivant de polynôme de Stieltjes et de 
Van-Vleck respectent l'équation différentielle fuchsienne : 


l=p j=p 


Zu racines de £ TI a =0 
He: l 1 
= y 
ye) Sul)=2-24 eil EOE -0 
l=p-1 E (z-a,) 


Polynômes de Heine-Steilties de degré 2 et Polynômes de Van-Vleck associés pour une «équation 
fuchsienne générale » 


Prenons maintenant le cas un peu plus compliqué de n=2 et p quelconque : 


SC Et né 
; | z ; 3 ; GI X=Y 
x, y racines de S,(z)=(z-x\z- y) Système d'équations fractionnaires i 
5 + A 0 
E yox 
Izp j= lz l= p» 
Posons D D AN Sek ] Il: oh SE et "tre Fab) 
m X74; x-a.) j=l I Xi age (x-a) 
j=1 ` l=1 i 
lži 
l=p_ j=p | l=p j=p 
P Jl Vi (x-a,) PO (x)= "7, (x-a,) 
"Je o Zu 
j=p j=p 
P (-2,) kal 
Sn o, 2 =0 Ge HN y=x+2À y-a=x-a,+2 À 
DIr Ah degt? Lal P (x) P a(x) 
j=l ! 
e NS) 4 E D =0 


Une mise en facteur commun de l'expression, un peu fastidieuse, conduit au numérateur à un 
polynôme en x de degré p(p-1), par exemple pour p=3, de degré 6 et pour p=2 de degré 12, comme 
suit que l'on va noter P,-1{X) : 


l=p y l=p y j=p ip rg, l=p j=p ip. rg, 
14 7 felt Als) =0S[[12+97 HS nii 2+} 7 ER 
Di G-a)xS Yj 12 E] 2+9 y A j=l i= =a; i= i 


ja X7 á; ii Xq 


P 
| ; ; ; : x-a p-l 
Comme BEE Sy A et Tr heh LE d Am IN Jl 
ii ZO Hi x-a; à X-a, [IG-a) Ai j| X—d; j=l 
` 
J=P P y l=p Jp i=p 2 PU) (x) l=p j=p PU) (x 
&|112+7, (x ob ` nii [| 2+7; (x ob i =06[[ 2+7; (x a) 5 v| [| 2+7; (x a) SS 
S SC A Gate [é=a)l "| à Tal 
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mal \ jel \ izj 


o ipfe «es 


mél \ jél ik) 


7 j=p p- i ` z (p-1\(p-2) 
(fie-e) D | 
ri GE nes Se | lee (ent 


Numérateur uf t uk ol (x 10) Se ii rel =0 


Numérateur — Polynôme de degré Dale p=3> fe p=4> x? 


La seconde racine y est toujours déterminée par l'expression : 


=p 
(x-a,) I=p j-p 


Fe ER avec P =n e-a) 


Il j=l 


Cette racine est toujours en correspondance unique avec l'une des racines x que l'on choisit pour le 
polynôme de Stieltjes. Il va de soi que l'on aurait pu commencer par déterminer la racine y qui suit 
exactement la même équation polynôme de degré p(p-1) puis la racine x en relation, comme suit : 


=l l=1 m=l j=l i=l 
+j mal Zei izj 


ET EEGEN 


Mais puisque x et y sont déterminés de la même façon, alors x et y sont deux racines distinctes du 
J=P j 
même polynôme de degré p(p-1), de telle manière que : Te-a) IA Ces 
y=x+2 et x=y+2 7 
(x) P, (y) 
deux expressions sont des involutions et définissent exactement un couple de deux racines (x,y) en 
correspondance l'une de l'autre. Il n'y a donc en réalité pas p(p-1) possibilités pour le polynôme de 


Stieltjes mais pn Uu fait qu'il faut choisir un couple de deux racines. 
2 


SE 


Pour ce qui est des polynômes de Van-Vleck d'après le théorème de Heine il y a au plus 


n=2 CMP? = Me 2 _ Px(p=1) formes possibles. Ce qui correspond exactement au nombre de choix 
nl{p-2)! 2 

de racines du polynôme précédent. Pour le calcul des racines tx du polynôme de Van-Vleck, cela fait 

intervenir les racines de la dérivée première du Polynôme de Stieltjes. Cette dernière est 

relativement facile à calculer à partir du choix des deux racines x et y de ce dernier. Notons là x’, 


alors: $()-(2-x{2-y)= 5,"(7)-2z GEHE "Let, Par ailleurs les p-2 racines tx du 


polynôme de  Van-Vleck respectent toutes l'équation fractionnaire suivante : 


£ 1 of? 2 ; , ; R 
ke12,..,p-2 + =0 + =0 
EE GEESS > EF Or les équations fractionnaires qui 


déterminent les racines x et y du polynôme de Stieltjes sont les suivantes : 


Donc Les deux racines choisies x et y du polynôme Pak) Sont bien solutions de l'équation 


e X i K l=p Ca 
fractionnaire des racines du Polynôme de Van-Vleck : NI , 2 -0. Désignons par ta, 
mit 2, -(x+y) E 


toutes les autres racines de cette équation fractionnaire liée au choix (x,y) sans être x et y. Alors on 
peut réécrire de manière formelle sans l'indice cette équation fractionnaire, en effectuant une 


dE i i 


mise en facteur commun, comme suit: Vi, 2 -0 j=l 2 


T t-a, 2t-(x+y) TK X2-(x+») 


j=l 


sachant que le polynôme au numérateur se factorise toujours par le produit des racines. 


Il suffit alors de choisir comme polynôme de Van-Vleck le produit des racines autres que x et y : 


Di | pour un polynôme de Stieltjes ` ,{,)_{,_;{,-,), comme suit : 


j=l 2 


l=p 


vc, Sala Jre guf es Joel -a (2+8)=0 614 Si =2+a+f 


Fr l 


Par ailleurs on sait que pour n=2, la contrainte sur æ ou 6 est la suivante : 


l=p Izp 
Terme z” »2+2x Ÿ y,+a B=0 et Dy=l+a+B=a=-2 ou B=-2 
l=1 l=1 


Donc l'équation différentielle fuchsienne est bien respectée par cette construction. Pour dénombrer 
le nombre de choix du polynôme de Van-Vleck, étant donnée que sa forme est directement 
construction à partir du choix (x,y), il y a bien p(p-1)/2 formes possibles. 


Un résultat partiel pour les polynômes de Heine-Steiltjes de degré3 et Polynômes de Van-Vleck 
associés pour une «équation fuchsienne générale » 


Soit les trois racines choisies pour définir le polynôme de Stieltjes : s(z)=(2-2X2-2,X2-2,) 


= 2 2 

D ' =0 
112474 un Zus 
l=p 

GM 2 2 _0 
12274 273 %77 
Si 2 2 

n =0 
Za Zar 735Z 


Ce système d'équation fractionnaire est loin d'être aussi simple à résoudre que celui pour les 
racines du polynôme de Stieltjes de degré. En cela le résultat ce petit paragraphe est donc partiel. 
Tout juste peut-on supposer que les trois racines sont solutions chacune d'un polynôme de degré 


p(p-iXp-2), et que chacune d'entre-elles est reliée univoquement avec les deux autres, ce qui 
2 


conduit à n'en retenir que le tiers soit p(p-1\p-2), ce qui est exactement les résultats de 
6 


dénombrement du théorème de Heine sur les polynôme de Van-Vleck. Car nous allons voir que le 
procédé de choix de ce dernier polynôme est similaire à celui réalisé dans le cas n=2. 


La dérivée première du polynôme de Stieltjes s'écrit : $:(2)=(2-2,{2-2,)+(2-2,X2-2,)+(2-2X2-2,) €t la 


dérivée seconde s'écrit :5"()=2-2, 42-24-27 +2-27,+2-2+2-2,-62-2(2+7,+-,) Les racines de la 


dérivée première sont celle du polynôme du second degré : 


2 
S,'(z)=3z 2z(z HZ, + zl H (z2, + 2,2, + zz) 
1 1 
ZE? 


S'()=36-7 2-7.) 


_ 22 +225 + 2223 


3 


z'+2!,= 3 (z, +2, +z) 


Les racines du polynôme de Van-Vleck respectent l'équation fractionnaire suivante : 


g Ir à 1 , 1 ` sl Vis 214-2727 
} =4 r Uc i J 7 0 
il h OZ IZ: md zi —-2; 
22! _ ZZ +223 +Z, 
= 61,-2{z +2,+z EZ: 3 
Di Vi k (z 2 :) =04 


lé 34? 21, ł t + t t 2 
De k k (z SCH z) ; (zz, E 22) rite (a+a+2) 


On voit que les trois racines choisies z:z, et z; du polynôme de Stieltjes sont bien solutions de 
l'équation fractionnaire des racines du polynôme de Van-Vleck (j'ôte l'indice des racines t,) : 


l=p l=p 
fes Vis = 6z-2(2+z,+z;) 0 Yi L2 2 zi=2 —2; 0 
SE } H 
Of Zi ZZZ -ZZ +22 A 274 (z DD z) 
ip 2 2 
en A i =0 


Comme l'expression est totalement symétrique en z:,z: et z3 les deux autres racines respectent 
également les équations fractionnaires de Stieltjes. L'équation fractionnaire des polynômes de 
Van-Vleck conduit par la mise en facteur commun à l'annulation d'un polynôme de degré p+1. 
Comme trois de ces racines sont zz: et z; il suffit de prendre les p-2 suivantes pour construire ce 
polynôme. Comme elles sont construites à partir des deux racines z;' et z,' notons ces p-2 racines, 
bis liée au choix (z:,z:z:) sans être z:,z, ou za, Alors on peut réécrire de manière formelle cette 


DNH 


équation fractionnaire, comme suit : 


l=p 
6+3 t t t Et, ; 
Yi 6t-2(z,+z, +z) 8 | Zo gd À 2) 1 sé) -0 
i 2 L L L L L j=p 
ia Milan le +22 +22) Gr 2t(2+2,+2,)+(27,+22,4 zz )x[ Ík-a,) 


j=p-2 
j= 


[I PR: ` Ni — lr z) 


sachant que le polynôme au numérateur se factorise toujours par un produit sur ces racines. 


Il suffit alors de choisir comme polynôme de Van-Vleck le produit des racines autres que z:z, et z3: 


j=p- 


Il ) pour un polynôme de Stieltjes : y(-)=( 


j=l 


z-z, z-z z-z) Comme suit : 


"cl ëch DEE CEA (2-2) ir ; E 
=p-2 fi LU l ds JA a EE Ss = = 
ra gf aT = (1) 6z ds es) 1: 22(2+2,+2)+(22,+22,+2,2,)+a p dë (z zz zl z) 0 
jH 1 
11 
Ep vr 
Ge kën Z-t, a) 
Comme G 2 -22(2+2,+2,)+(22,+22 + 2,2,)x "i 6z 2(2+2,+2)= S = Gaza 
Z-A 


BE E TE ET DEE DEEN 


[IE-a) Lal 


ja II 


SeS y ta p=6+3(1ta+ fra B=06(+a)5+8)-0e Fy, =l+a+p 


I l=1 


Par ailleurs on sait que pour n=3, la contrainte sur o ou 6 est la suivante : œ =-3 ou p--3. Donc 


l'équation différentielle fuchsienne est bien respectée par cette construction. Pour dénombrer le 
nombre de choix du polynôme de Van-Vleck, étant donnée que sa forme est directement construite 
à partir (z:22,z3), c'est le même dénombrement que pour le choix des trois racines. 


Par extension on procède de la même manière pour un polynôme de Stieltjes de degré n, même si 
les expressions sont plus complexes. Les expressions des dérivées première et seconde sont les 
suivantes : 


i=n i=n i=n-l i=n i=n-1 i=n-1 


s,G)=TIt z) E S,"(:) (t-2',)=¥ (n-i) = Hi Us, Al =nx KS TT ` ei, produits symétriques 
i=0 


i=l i=0 i=l i=0 i=0 


où set s' sont les produits symétriques suivants : 


l=n i,j=n i,j=n ol l=n 
s =l H Z sS, = Da s; = dre dé ST ge, da E Zi 
1 ij=Li> ji i,j,k=1,i> j>k iset EE l=1 
l=n-1 i,j=n i,j=n DË l=n-1 
fa RS H Lë ta Š 1 St = ' ' F a ' 
DER SAS Zi E KEE s, = KEE Sé KIT KE Za Zi 
1 i, j=li>j i,j,k=1,i>j>k DEE Ed l=1 
i+ 1) A n-i 
' EE ( n-l-i ' 
H n-l-i Ès H 
n n 


i=n-l 


La dérivée seconde du polynôme de Stieltjes s'écrit :5"()- ing.) et nous avons: 
i=0 


SR Si bien que l'équation fractionnaire des racines du polynôme de Van-Vleck peut 
Sz) 2-7 
s'écrire : 5 n ell ua 1, lu. Cela conduit à une équation polynomiale de degré 


df t-z med ët) | AU 
p+n-2, dont il faut exclure les n racines solutions correspondant à un choix de n racines du 
polynôme de Stieltjes. En effet ces dernières sont bien solutions de l'équation fractionnaire puisque 
l'on retombe sur le système d'équations fractionnaire des racines du polynôme de Stieljes. 


t= A > d Vi f AO = g Vi Ste) 
l=1 t-a; Si l=1 Zi =a Ss 
" l=n i= l=n 
Or S, (z) 2 2 > £ Yi > 2 d 0 
St Zä 11 Zi dj i74 


Aussi il convient donc bien de choisir pour racines du polynôme Van-Vleck les p-2 racines qui ne 
coïncident pas avec n racines z;. Notons ces racines comme suit ` et le polynôme de Van-Vleck 


GE j=p-2 ÿ k 3 R 5 
ainsi y(z)=a P Me-t) . L'équation fractionnaire de Van-Vleck par une mise en facteur commun 
j=l 


conduit à l'équation polynomiale formelle suivante : 


l=p l=p y l=p l=p i=n j=p-2 l=p j=p-2 
SEA sO e-a) [$n -dhie ffe- Kl Bä aleet D 
l=1 l1 t74; 121 Ial gy ja m E sa Bogis 


On peut à profit utiliser cette factorisation dans l'équation différentielle générale qui se trouve bien 
vérifier : 


j=p-2 
sl: S, (z) l=p y E 
E END L jxS '(z)+a p xS,(z)=0 
V(z)=a DUT d 5 ge Ss fi (z-a,) S 
j=l 1 
IA 
l=p l=p l=p l=p j=p-2 
come Se} -as offe- 5- mn satt, 
l=1 l=1 l1 Z7 IA jal ES 
En 


usA y +n(n-1)+a p =nx(1+a+p)+n(n-1)+a B=(n+a]n+8)=0 


I 


Construction des développements de Frôbenius autour d'un point singulier z=z 
Recherchons la solution de l'équation différentielle suivante : 

y"(z)+ P(z)y'(z)+ Q(z)x(z)=0 
autour de la singularité z=z, constatée sur la fonction P(z) par un développement de la forme : 


j=+0 


Z@=(-2) 3. CG-2) 
j=0 
Dans la nomenclature de Bôcher, l'équation prend la forme suivante : 
Z"(z)+ P(z)Z'(z)+Q(2)z(z)=0  (z-z)P(z)=% 4-2) (z-z Q(2)=2, B-z) 


Introduisons l'opérateur différentiel fonctionnel suivant: 


L=(-2Ÿ S 2 +e- seat +e fat 


P(z)={(z—2,)P(z 
Se ` steet Zeckt 
Olz)=(z-27 Q(z) 
L'équation différentiel peut alors se réécrire sous la forme : L(Z(z)) 50; 
En prenant une solution sous la forme de la série proposée, l'effet de l'opérateur différentiel L est le 
suivant : 


zen= È C e-2)" > Lars ECS) È Eeler 


Autour de la singularité z=z,, on développe les fonctions P(x) et Q(x) : 


i=+0 


fe-zoor+i)=fr+ il +[r+ dn e (+ 5Xz-2) 
ee zl J-E cn (r+ jX z-z ea 


En réarrangeant la sommation de telle sorte que l'on rassemble les mêmes puissances de z-z0, il 
vient : 


L(Z(z,r))= De A) "AC f(r+v)+ Cfifr+v-1)+..+C fr) 


SLGGr)= Alsace Cf +v-1)+ le zu 


v=0 


Si la série est solution de l'équation différentielle, alors chaque terme de même puissance doit 
s'annuler ce qui conduit à résoudre un système d'équation linéaire sur les coefficients. 


yv Cfir+v)+C, fifriv-1)+.+Cf(r)=0 
Chi(r)=0 

Cffr+1)+C, fir)=0 
Cffr+2)+Cffr+1)+C f(r)=0 


Cfr+v)+C, fitr+v -1)+..+C,f(r)=0 


C'est la forme très générale de solution en dehors du cas de l'équation de propagation de Bessel. 
On peut également comprendre ce système linéaire comme un système de récurrence sur les 
coefficients, comme nous venons de le faire dans le point précédent sur l'exemple de la fonction 
d'onde de Bessel. Dans le cas général, les coefficients se déduisent par l'expression suivante : 


r tq hr)=r+(4-1r +B =0>re{fnn} n>n, f(r)=r4,+8, 


al) -c Yal) 


REN A E "rett 
(OH Da ou TE 
he) Ae+1)  ffr+2) 0 0 
ul Ar) är gead 0 
e E f,atr+1) f,a(r+2) fal +3) hr+j-1) 
Ziel Zei f+) ffr+3) - fAe+j-1)]| 


Comme le déterminant en question est identique si l'on range ligne et colonne en ordre inversé, 
alors le déterminant peut aussi s'écrire : 


ët he+j-2) ffr+i-3) ~ faei) Zill 
HU EE) LE) à oea 

8 (-1¥4 (r) se 0 f(r+ j-2) fitr+ j+3) e, Ge 
SE E Pa Air) . 0 p vi A 2 
SE se 1 Sc? fr +2) fi(r+1) hlr) 
0 0 0 fAfr+1) Ae] 


A(r) = F, (r) Notation de E.L.Ince 


C'est exactement la notation et l'ordre du déterminant utilisée par E.L.Ince (page 397). 


Tous ces coefficients et tous les déterminants doivent se comprendre comme étant des fonctions du 
paramètre 6. Par exemple : 


z S TAE) he) TATAG. ` 
EE EE A 


Si nous reportons les valeurs des racines de l'équation indicielle rı et rə dans les formules, il vient 
pour le dénominateur des coefficients : 


Girls +(4-1r +B, =(r-r r-r) En E E E Den pour Dénominateur 


(- 1)’ A; (r) 


z , ! T-r+j+1) Tfr-r, + j+1 
= Ci(r)=C F et Den(r)=] [(r+i-r) (r+i-r,)= toits 3 Kë OI 
(r+i-n\r+i-r) gë E l 2 
i=1 

pen) Me » 8 
n>r,n=n+k> | : ou bien Den: )=] [ix] [C+4) Den(r,)=] [ix] A 

Del )= DCL E +J +1) i=l i=l il D 

| T(1-k) 


En d'autres termes si les deux racines ne sont pas distantes d'un entier, et k restant un réel alors les 
dénominateurs sont parfaitement définis. Et l'on peut alors écrire : 


Li (—1Yr(k+1)A (7) bit (-1)T(1-K)A (7) 


JT (k+j+1) °  JIE(j+1-k) 


Les cas de figure se compliquent lorsque les deux racines sont distantes d'un entier k >0 car la 
fonction Gamma présente un pôle, ou bien que les deux racines 68; et 6- sont égales et dans ce cas 
par construction les deux solutions coïncident, il faut donc élaborer une seconde solution 
indépendante de la première autrement, et pour cela opérer une légère modification de la 
méthode. 


Modification de la méthode de Frobënius pour obtenir les autres solutions dans les cas des racines 
indicielles multiples ou séparées par un entier 


Voir l'ouvrage de E.L.Ince « Ordinary Differential Equation », page 397 et suivantes (point 16.12). 


La série introduite ici a une forme similaire mais le paramètre indiciel n'est pas racine de l'équation 
indicielle. De ce fait le développement devient une fonction dépendante de deux arguments r le 
paramètre indiciel qui devient formellement continu (discret lorsqu'il est racine de l'équation 
indicielle) et z un point du plan complexe. 


Z(zr)=(2-2) F C, (rXz-z,) C(r)=C, quelconque et Girls 
j=0 


Toutefois les coefficients de la série sont solutions du système linéaire qui suit : 


Gr +1)+Cs(r)=0 
Cr) hlr +2)+C,(r) (+1)}+Cf(r)=0 


lt suis Ctr)ffr+v-1)+..+C(r)f..(r+1)+C, f,(r)=0 


Alors l'opérateur différentiel appliqué sur ce développement en série a une valeur très simple 
puisque tous les termes de puissance positive s'annulent (équations du système linéaire). Il ne reste 
que le premier terme : 


Ce résultat n'a l'air de rien (E.L.Ince, « Ordinary Differential Equation, page 397), mais il est en fait 
très utile pour la construction des secondes solutions. On retrouve également la construction des 
premières solutions si 6 est racine de l'équation indicielle : fo(r) =0. 


Cas des racines entières distinctes ou de racines réelles distinctes dont l'écart k est un entier 


Donc lorsque k est un entier le terme r( j+l- k) présente un pôle, on peut formellement séparer 


les coefficients du développement comme suit : 


e j<k Ur att Ara a 


rei a ai nb. — 


On a donc isolé le terme de divergence des coefficients d'indice supérieur ou égal à k. En reprenant 


À A. A 
mot pour mot le raisonnement de Moon et Spencer, on peut supposer que le terme Ain soit 
aie r>r 
à tout hasard fini, si le déterminant s'annule sur les deux racines indicielles, mais comme il ne faut 
pas trop compter là-dessus, on préfère introduire pour ce groupe de coefficients la grandeur 
suivante D;(r) qui neutralise l'annulation : 


7 +1- H 1} j Ee) 
E Uu ne k) Ge BK 


i=l 


D(r)=(r-7,)C;(r) pour j2k D,(r)= 


Pour j>k D (r,)= Lim, „(r-r )C,(r)= 


k Î 
olli, -r ci (r-n +i) Dig WIIG- 


i=l i=k+1 i=k+1 
GOUT CA () 


GG -E) 


k k 
Pour j=k Dj(r,)=Lim,.,.(r-r,)C/(r) CU CA) ` KU Giel _ C, A(n) 


> D,(r,)= j>k 


On complète la définition des coefficients D;(r) pour les indices inférieurs à k. Par construction 
comme les coefficients C{r)ne présentent aucune singularité, les coefficients D;(r,) sont donc tous 
nuls : 

Dirk G—nkKicl- kk Lim... nich 0 pour j<k 


Ces coefficients servent dans la construction d'une seconde solution, dans le cas des racines 
entières distinctes, sous la forme : 


= Ž,(r,z)= let +20 earth D Ditz" 


d 


Selon la théorie de Frübenius, on parvient à construire une seconde solution indépendante de la 
première, par un processus limite de dérivation sur la seconde racine indicielle, soit : 


> Z,(z)= Lim, An seconde solution 
r 


IST 


Par ailleurs dans le cas de racines indicielles séparées par un nombre entier, il existe une relation 
entre les déterminants utilisés pour calculer les coefficients Ci(r) et Dir On peut tout d'abord 
montrer que le rapport des coefficients D;(r) tend à la limite vers une expression définie. Ici on 
prend la précaution du passage à la limite car les coefficients Ci(r) sont censés diverger pour r= r2. 


Cl bel Cr) 17 KIA (7) 


Jim, T T TE 
ot) Dh) C, Jj +k} 


Pour établir cette relation il suffit de considérer le système d'équations linéaires permettant 
d'établir les coefficients D,.(r }par passage à la limite. Commençons par les coefficients C;.{r) à 
partir du coefficient C,(r ) et appliquons la définition des coefficients D;.(r ), il vient : 


Giele A r+k+1)+C, (r Mie+k)+.+Gfr)ffr+1)}+ Oh fial) 


de E r) Atk j-1)+.+ COS elt Cf, del 
Si r=rn r,+k=r et Lim, (r- r, )C KG )= Dir 


Dalali +1)+D, GE e - 1)+.. + Dill -1)+ O(n -r)i -k)=0 
D, lte ep. MA ++ Dh Ae + Gr) 20 
Comme Dfr)=D,(r)=..=0D,.(r)=0 
D) +1)+D,G)A(G)=0 
_ Peal) or +2)+ Deal )Al +0) + D (5) (5)= 0 
D, elt +) Da elle +-+ D) Bl) 
DJ E EE E le 
goe ép 
D.A D q Dou £ En Ce 
= D,(r,) Lin 2) D,(r,) filr 1) Gin 0 
D Jl Lu RA GE 4 Dalh) Mers 
D, (r) Lin j) Dik At J 1) Wi DE” i(n- 1) CAE 


Ce dernier jeu d'équation linéaire est exactement celui formé par les rapports Ci(r:)/Co à savoir : 


"luet ig 

Gel cp 21 El ob Ae (= 0 

C, oI C, 11 2V1 

ln Ell pr js El. adese) 


C, C, C, 


Autrement dit les rapports étant calculés par le même système de récurrence (ou système linéaire), 


ayant le même point de départ, ils sont égaux. Soit ol) = Pers (z) Vj>0- 
C, D, (r,) 


Ici je reprend le raisonnement de l'ouvrage de P.Moon et E.Spencer « Field Theory Handbook » à 
ceci près que pour ôter toute confusion, il faut utiliser les nouveaux coefficients D;(8) en lieu et 
place des anciens coefficients C(8), ce qui n'est pas indiqué dans le livre et cela rend impropre 
quelques formules reportées pour la démonstration même si le résultat final est juste. 


Et par ailleurs la relation sur les Dir peut également être établie comme suit : 


GA) Au. e p T A 


D (r,)= — 
1) AG -DG-#) DI ai e EN 
(-1V"A a(n) 
— Din) _ jI(j+k} = _1} KA a(n) SEN Call Dall + y KIA aln) 
D,(r) _ GA) JG+4)JA,() té C(r) D,(r,) JIJ +kXA, (r) 
k! 
Les deux expressions comparées, il vient : 
, Dain . Ha, Alr) 
E E GD j+k V2 = À. Dirk V2 
ge tea) AAO fr) 


Soit la formule importante sur les déterminants : 
All A,(r)A,() 


La dérivation de tous les coefficients du développement lorsqu'ils ne présentent pas d'annulation 
du dénominateur est la suivante : 


| (- 1)'A (r) 


(r+i-rXr+i-r) 


Comme C, (r) =C, 


1 i=l 


Di ac (r) EU GT DE ail Ee Seil 
) 


"mn Gr JIT(k+j+1) | Or 


GC) GI) [2k Al | 


P i<k Li = 
our j< Im. Gr Gan 


Pour les coefficients supérieurs ou égaux à k, on dérive les coefficients D introduits précédemment, 
comme suit : 


aD, (r) eG ôr 


D ` Oil uff a D = | 


i=1 il EA Sekt, te Ee Et 


YO 
D. GE Eil 
D ee ` Gr: SL. A H 
NTI 
Or Ar 
k 


vk | j k-1 j- 
sims Og CD CH Ain Ze Së d 
I=1 


"og "ott 


Nous avons vu comment dans un processus limite et dans le cas d'un écart entier entre les deux 
racines de l'équation indicielle, construire une seconde solution 


Nous donc maintenant considérer les trois cas de figures. 


Cas des deux racines réelles distinctes d'écart k non entier 


Les deux solutions de première et deuxième espèces se déduisent des coefficients calculés 
auparavant qui n'ont aucune singularité au dénominateur : 


en cc) pee EIS 


TE + j +1) jT(j+1-E#) 
At) ` de) 0 0 o | 
AC) Oezl pesa 0 0 
TE 
Ía r) far +1) fia(r+2) fi-a +3) he+j-1) 
L f(r f.(r+1) f.(r+2) f,s(r+3) filr+j-1)| 


Cas des deux racines identiques (réelles ou entières) 


La solution de première espèce est identique : 


k=r-r=0 ==) CC Ain) 
! (j!) 
ft) ` Gs 0 0 0 | 
ft) ft) ffr+2) 0 0 
a Isi. ft) A+) E f(r +3) 
La) fatr +) f(r +2) far +3) hr+j-1) 
r Fae lee Eh, e=] 


Pour la seconde solution on introduit une fonction de la forme suivante, en considérant la racine 6o 
comme une variable , soit : 


En introduisant l'opérateur différentiel fonctionnel L et en passant à la limite sur la dérivée du 
paramètre de racine indicielle, il vient : 


2 


1-2) +22) IE ka Fast 


z z 
Si rn > f{r)=(r-1) 
E (r, z))= C r-n (2-2) = Cfr-r,ÿ e' #6) 


E ô tene d aa 


ôr Or Or 


SE e 
r 


- JS E — Ae = 0 


De plus : © 
ð 


Z,(r,z)n'est pas une solution de l'équation différentielle 
r 


Mais en dérivant par rapport à la variable racine indicielle et en passant à la limite, la fonction 
obtenue satisfait à l'équation différentielle, il est donc théoriquement possible de construire une 
seconde solution, comme suit : 


. ô z A EES 
Z,(z)= Lim, EL (r,z) constitue une seconde solution indépendante de Z,(z) 


Or 


Il vient : 


r Je r+j 
Z,(r,z)= Cz- 27 + C,(r\z-2) f 


j=l 


Z, (r, z) ES Ce Log(z-20) +e ee Ç; (rXz Lg y 


—Z ( ) C Log(z z aJe Ek: w + Log(z- Zz aJe SR ka c Xz- zl EE 
j=l j=! à 
= C,Log(z- z (z- 2) + Log( Gs z) IS ci Xz- Tay EC 
j=l j=l d 
Wal) . oC o YT aA) mi 
! = = JA. = 
Ci) Co (GY ES C, (GY CH An, H I 
= (17A, 
C,Log(z— zz- z) +C Logis a): l ) am DH A 


J 
1 2 
=> Lim Le (r,z)= a 14 


ee eyi a oaa a S 
SE (JY f ôr Al ` 
ol 


Ze (-1)'A, 
Comme E ERS le 2 D F SEAN 


i=» OO AT Fo +) | l=j 
> Lim, ENE E E C, ES ( 1) e Ri ak 2A (r,) SE 
° Or = (j) d j 


Soit pour seconde solution 


(Gef AM ro + OA. l=j 
Sien EE e l (n) 24 (82) 
KR Or j=l (Lal Or lz J 


Il peut arriver que cette construction ne coïncide pas avec des secondes solutions usuelle dans la 
littérature scientifique. Dans ce cas, il suffit de rajouter un terme proportionnel à la première 
solution pour faire coïncider les deux solutions, comme suit : 


EE | Fi EE +j ðA de ) tj 1 
Z,(r,z)= Aë zl Zi(n,z)Log(z 2): C >, l ) l 5 2) l Se ELO 
Ja (j!) ôr mi j 
La seconde solution se calcule donc sous la forme suivante à partir de la première solution avec un 
terme additionnel, avec éventuellement adjonction d'un terme proportionnel à la première solution 
pour recouper les définitions usuelles des secondes solutions : 


Z()= Loglz -220+ Gl-2) È z CH 2) E (a d e ` 


Cas des deux racines entières distinctes ou de deux racines réelles présentant un écart entier 


La solution de première espèce est identique : 


k=r-r#0 keN r#r=C{fr)=0C, 


! ji(k+ j} 
Ae) ftr+1) 0 0 o | 
ft) ft) f(r+2) 0 0 
ae O ëch SH) ar 
AN D) at) 0) rar AA) 
FO) HAN Le) For) sr Ar EE 


Comme il est indiqué dans le processus de construction de Frübenius, on introduit une fonction de 
la forme : 


j=1 j=k 
Comme D,(r)=(r-r7,)C;(r) pour jzk>Ž,,z)= Y be )C (rez E? 
j=0 


où les coefficients C6) sont solutions des équations linaires : 


Ci(r) Je +1)+Cfi(r)= 0 Cr) o(r+2)+Ci(r) (r+1)+Gf(r)= 0 
Lu (r+k-1)+C, (r)i +k-2)+...+C(r) allt Chr 


Get C, (r) (r+v -1)+...+ C (r)f, a(r +1)+ C,f,(r)= 0 v2k 


En conséquence l'opérateur différentiel fonctionnel L a une valeur différente de celle établit en 
début d'exposé sur la théorie de Frübenius. Cela est dû à la construction de la série utilisant les 
coefficients D;(r) au lieu des coefficients Ci(r) à partir de l'indice k. Toutefois il reste toujours 
uniquement le premier terme de la série, non nul : 


dëch aleet: E Ce-a)” E Deea" 


n Zens al" Cea) Cae). Geh 


v=k 


= L(Z(z,r))= Cz- 2 (r-r) sel Cr -ny r-7Xz-2) 


Si maintenant on dérive la fonction par rapport au paramètre indiciel, et que l'on passe à la limite : 


[EEA] -Ecne 
Or Or 


=2C, (r -r, Xr -7 k Log(z-20) +C, (r -r e ef Log(z-20) +C; (r -r, P (r -r )Log(z A k Log(z-20) 


Comme dans le cas précédent, la dérivée par rapport au paramètre de racine indlicielle, satisfait à 
l'équation différentielle. Il est donc théoriquement possible de construire une seconde solution, 
comme suit : 


Z, (z) = Lim 2 Z, (r, z) constitue une seconde solution indépendante de Z, (z) 


FT ôr 


Plus précisément la dérivée de la fonction introduite s'écrit : 


C, (z — Zu Y Log(z SÉ )+ 


` C, SCH el dE (r) 


aZ (r,z)_] = +5 Pic peste -a)le-2)" 


Or | j SS 
Jato »+5 | 0D; (r) 
+ (z-z) —. t D,(r)Loglz-2) 
Jk 
CJe-z Y Log(z SÉ )+ 
O Li ror, 7 GH ôC, +j =0 
Se ee + i E FC, (r)Log(z- z, le -7 
Ja 
7 j=k-1 | j=+0 | D. 
= kinn PEDE eier Ein o (note) 
j=0 j=k 
E EH ôD, (r) (- EE ôA (r) iq HI Di 
O D, = C A , = C A d __— Ee = 
Fe) (att NZ e (att 001 o MEA 
Gärt 
leche, T A (r) j<k 
Ciz-aAf © (k-j-1) j 
Ve D F A (r (z-z) + 
C (z -zZ T de hä (z -z y 
Z = E 07 La 
= SEN + (Gu og(z gd GE) (2) 
Sab tre E E sa 
St Ster et 


(Gu & oU ôr 


Toujours mot pour mot de l'ouvrage de Spencer, la deuxième sommation peut s'écrire par 
translation de l'indice j, partant de 0. Pour cela on utilise la relation suivante sur les déterminants : 


A jak (n ) = A, (r Ja, (n ) 


Il vient donc : 


Clz-2) © (k-j-1} j 
Y 2, F A (r z-z) + 


Alt toute) (EEE, Za aal 
Glz- CNY all län) LI 1 SI 
“He e "ER | viel Sr d 
Comme EE 


bastele LES EE ae) 


alea P CNE aY D Al 


(Gu &  jI(j-k} ôr 


Il se trouve que l'on peut également décaler la dernière sommation de la même manière : 


A) z (hrog(e 2) ER SEE (je) 


SEAN ELE 


(Gu & oz ôr 


> Z,(z)= 


+ z(e) f DEA SEE Lat 


Gle- n) S (1) (-2,) Per ke ul e ZO 


(GA & jj+k) ôr 


Outre la forme précédente que l'on peut également utiliser, la seconde solution se calcule sous la 
forme suivante à partir de la première solution avec un terme additionnel : 


Log(z Zi A (z)+ 


A(n) 
DIE 1 


C lz-2) JS C 
(Gm 2 j! 


j=+e (- It 8 H ôA (r) à Izik] ii] ii 
z S AE pi SS 
(GAN &  (J-k)j! ôr e mi ml ml 


Z,(z)= + Z 2) + 
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L'approche de E.L.Ince : construction des secondes solutions dans le cas d'un écart entier ou nul 
également par passage à la limite 


Ici je vais tenter de reprendre les calculs réalisés par E.L.Ince, dans le contexte des seules équations 
du second degré. Ce dernier les réalise plus généralement pour une équation différentielle 
ordinaire de degré n quelconque. Il s'agit du point 16.3 page 400 de l'ouvrage « Ordinary 
Differential Equation ». Et je me place délibérément dans le cas de deux racines r1 et r2 dont 
l'écart est un entier. En outre r1>r2 et le le point singulier est z0 au lieu de 0. 


Selon les notations employées par E.L.Ince : 
Z(z, r) = ` C, Liz SEN E" Cy (r) =C, quelconque et hlr) +0 
j=0 


On a déjà vu que l'opérateur différentiel appliqué sur ce développement en série a une valeur très 
$ D A 
simple : LZE, r)= Glz -2Y EE 
Z Z 


Les coefficients sont construits formellement à l'aide de déterminants : 


fGr+i-1) Str ffr+i-3) — Saet) Zill 
y , , fhr+j-1) GE EE a f..(r+1) far) 
Chic, SA F (r)= 0 Le A+ j+3 cé d 
i d Së t sc W X ffr+1) flr) 
LIot+) R i a AD a e 
|0 0 0 hlr+1) flr) |] 


Dans le cas où les deux racines r1,r2 sont soit égales soit d'un écart égal à un entier. On a vu que 
les dénominateurs s'annulaient et que dans ce cas les coefficients devenaient infinis. 


Cas de deux racines r1,r2 d'écart entier k >0 


Pour éviter cela E.L.Ince propose de multiplier le coefficient CO par une constante qui va compenser 
l'annulation des dénominateurs. Si nous posons k=r1-r2, alors 


k=r-n CG—Cfr+1)f(r+2).f(r+#)=Cf(r) 


La solution est alors construite comme suit : 


ce-a PEO äi SO GE eNe) 


H 


IC +5) 459 (uc 


On note g, (r) =C, (r)f(r) —> Z(z,r) = 2 Z, LL 2, y” 


On s'assure maintenant que le développement est défini pour r=r1 et qu'il n'y a plus de divergence 
sur les coefficients du développement lorsque r=r2 : 


eue) Ï 


j=+0 


Du r+i) )=>Z(z, r)= `X Cl UN F z-z)” 
ge [LA +5) 

> Z(z,r)= > Co [LAC + F, Xe -2)" + > C, CRU 2)" 

fe D Toi 


i=k+1 


i S P E -1řF (r 
aen TIn) Zn Be, Iae- Se, EOG- 


i=j j=0 j=k+1 | 
i=j+1 i=k+1 j i=j+1 j=k+ l [hli +i) 


k j ag a —1)F (r 7 M 
[LAC +i)=0 [LAC +540 7Z(27)- Ci C ) Je (z Ed 


i=j+l i=k+1 j=k+1 II (r. +i) 
0 (2 


L'opérateur différentiel L appliqué sur ce développement donne : 


(Lëtz cl Cl- AOE- Ar +). fr +2) 
Notons F(r)= Hrltde (tes Kl LZ(r)= Ce af Eil 


Notons que Allen ke donc le facteur F{(r) est de degré 1 en 
en ee ae ms on 


{r-r1) et de degré 2 en (r-r2). Nulle part ailleurs dans le facteur F(r) n'apparaît autrement les deux 
facteur (r-r1) et (r-r2). De même EE et le facteur f(r) 


fhlr+k)=(r+k-nr+k-r,)=(r-r,\r+r, —2r,) 
de degré 1 en (r-r2). 


Par ailleurs l'opérateur différentiel L commute à avec l'opérateur différentiel sur le paramètre 
indiciel r. 


s=0 pour rn 


et conséquemment Cup =0 u=1,2 | est une solution de 
Ta 


r=r, 


s=0,1 pour r=r, 


l'équation différentielle de départ. 


j+ 
Zo) 


Si l'on écrit le développement Z sous la forme suivante : 


J=+0 Ï=+0 4 +00 


Z(zr)=(z-2) 2, C, (r Yf(rXz-2,) =(2-2,) 2 g,(rXz-2) — g” Telz) el Xz-2,) 
+ À Hana) e + Lonte-aXe-a)  50%e-n) 
AE w(z,r)+ Log(z-z,}m(z,r) avec w(z;r)=Z(zr) et w(z,r (z- z) D (rXz- zF 


Or 


Considérons la première racine r1, alors il n'y a pas de dérivation nécessaire et nous avons la 
première solution avec des coefficients toujours finis (ou nuls éventuellement selon la valeur de 


G(r1) ): 


ZGr=wmer)=G-2) A. gn- ` gisl-Gselrisl fini 


J= 


© 


Considérons maintenant la deuxième racine r2, alors une dérivation est nécessaire et nous 
obtenons la deuxième solution indépendante de la première : 


ES 
Z, (z,r,)= zer) = w(z,r,)+Log(z-z,)m(z,r) 


Récurrence des coefficients g;et H 
La récurrence des coefficients C;étant la suivante : 


vv Cfr+v}+C,  fifr+v-1)+..+C f(r)= 


Alors la récurrence des coefficients g; suit une récurrence identique : 


>g (AC +5)+8 CAC +5-1)+.+ gelt 


Sr) C > H()f (+ j)+H,,6)f(r+j-1)+..+H,f (r)= 


_ at) Gk) CF) 
go(r) Co E 


Cas d'une racine multiple d'ordre 2 : 


Dans ce cas les deux solutions construites sont strictement identiques : 


Sich Se, EU y 
JH (ici 


r 2 
Comme ZE) =C al Al avee Zellen et que l'opérateur L commute toujours 
avec l'opérateur 2 il vient : 
la 


s=0,1 pour r=ņ =r; 


JH 


= Cfl -7} Loglz-2,X2-2) kä (r-r) 


= Cl(r-5Xe-2) 


BE 


donc là encore la dérivée formelle PC n) est une solution de l'équation différentielle. 
Or" E 


La première solution s'écrit : 
Ze rj=men) Gales Se / C (i z-z C (n) fini 


La deuxième solution indépendante de la première s'obtient par une première dérivation sur le 
paramètre indiciel comme suit :: 


Ion Los 


E Chefs seat È CXe-2) 


ran j=0 j=0 


Points singuliers apparents et absence du terme logarithmique 


E.L.Ince déduit de ces calculs formels (voir paragraphe 16.33 page 404 et suivantes de e Ordinary 
Differential Equation »), la condition pour que les termes logarithmiques n'apparaissent pas. 


Tout d'abord pour que toutes les solutions indépendantes développées autour du point singulier, il 
est nécessaire que les deux racines de l'équation indicielle soit distinctes. Le terme logarithmique 


n'apparaît que si l'écart des racines est un entier. 


On peut donc arranger les deux racines r1,r2 avec r1>r2 et r1-r2 € N. Reprenons les raisonnements 
employés par E.L.Ince qui sont semblent-ils repris des travaux initiaux de Fröbenius. 


La solution liée à la racine r2 se construit comme suit : 


Ô í j=+%0 , Ve" 
SE =(2-2,) £ g; (r) E (z-z) +(z Zo)  Log( Z— Ze 13 g(r (z- 2) 
Tt j=0 S j=0 
La condition pour laquelle il n'y a pas de terme logarithmique et alors: vi g(r) j 


conséquemment étant donné que les coefficients gj sont essentiellement des coefficients 
polynomiaux en r, alors eil œ(r-r,)P(r). On sait que le facteur flr) est proportionnel 
également à r-r2 donc gO(r)=COf(r) est également proportionnel à r-r2. Il vient donc que le facteur 


H (7, )= gi) doit rester fini (ou éventuellement nul). Prenons maintenant un des résultats de la 
J 
DU 


récurrence vérifiée par les coefficients H; : 


E 


J H()f(r+j)+H,,6)f(r+j-1)+.+Hf(r)= 
[TAC +5) 


Donc par cette récurrence si les coefficients H,(,)H,(,)...,H,,() sont finis alors le coefficient 


j- 
(r hlr, +j) est également fini. Donc H (r,) est fini si et seulement si fhlr +j)#0, ce qui 
J J 

n'arrive que si j=r1-r2. Donc Hir est fini pour tous les coefficients excepté pour le coefficient j=r1- 


r2. Et la condition nécessaire et suffisante pour laquelle ce coefficient s'annule est manifestement 
que le déterminant TIET F,(r,)=0 s'annule. Donc la seule condition pour laquelle le terme 


logarithmique n'apparaisse pas pour la seconde solution est que le déterminant suivant soit nul : 


Aer DH) as BA) CN 
ft -1) Ge GE ie E fe (6) 
SES 0 fhi- 2) fi 3) Sé éi 8 
Lech 0 er," HE) hi 
S — ….  … ris Aln) 
0 0 … 0 ff +1) vil | 


Exemple illustratif de l'absence de terme logarithmique sur une équation différentielle du second 
degré 


Voir le point 16.401 de l'ouvrage e Ordinary Differential Equation » page 407, E.L.Ince. 


Soit l'équation différentielle d'ordre 2 ‘2 y"(z)-(42+2 2° )y(2)+(4-x z)y(z)=0: ZcH est un point 


singulier. 


Suivant les notations on a 


4 4 
net 00)-4-5> m= A 
L'équation indicielle est donc la suivante : p =-4 q,=4>r-5r+4=0 (r-4Xr-1)=0- Les 
fonctions fi sont définies comme suit : 

= : fr =(r-1X{r-4)-2(4 r+K) 


f(er)=r(r-1)+r zP(z)+ 70 ()- ZA = \ftr)=t-1Xr-4) Aer) 
äis 


zP(z)=-(4+2z) z Q(z)=4-xKz )=0 i>l 
La relation de récurrence pour les coefficients du développement est donc la suivante : 


C hlr+v)+C n fr+v-1)=0 6 C(r+v-1)r+v-4)=(2 (r+v-1)+K)C,, 


C'est une relation de récurrence à deux termes. Les deux racines indicielles sont r1=4 et r2=1, soit 
un écart r2-r1=3 qui est évidemment un entier. 


La première solution pour la racine r1=4, s'écrit sous la forme : 


y(z) = cz“ (1 +yz +y,z? + a) 


Y,  ZDGrvbs be 


e vG + v) 


Les déterminants s'écrivent assez simplement comme produit de la diagonale : 


[flr+j-1) 0 0 [ 
AG+i-D f{r+j-2) 0 7 
F (r)= 0 eo on i S ` 
ei e Li … fi{r+1) 
| 0 0 ep 0 fifr+1) fir) 


> Eist AA (+). feu 5-2) (r+ 7-1 


La condition pour laquelle le terme logarithmique disparaît est donc la suivante : 


F,,(5)=0=F()= All) 
e Unie, +1)+k)(4 (r +2)+K«)=0 


K =—À 
&(2+k/22+rk)31+k)-0e {x =-21 
K=-3À 


Dans ces conditions la solution sur l'indice r2=1, se développe ainsi : 


En déterminant la récurrence générale sur les coefficients g'i(r), on obtient les développements qui 
s'avèrent finis : 


PA= ==) 
ar)=Cftr+0)f(r+2)f(r+3)=C(r-3Xr-2X7- Dr +1Yr+2)= g,'(1)=12C, 


Posons C,= = > g'(1)=1 


"(r)= '(r Alexis, CC 
(o) ( etj- 8; (r) 8 (r) FET) Zal a(r) 
8; £a fr + j) SE d É Les 
dr Ale: 


g(r)= Cr -2)r-1)r +1) +24 r+Kk)= a (1)= 0 
g(r)=Cfr-1{r+2)2r+KkX2 (r+1)+x)= g,(1)=0 
K=- > g,(1)=0 
gr)= GG r+kXa (r +1)+e)4 (r +2)+ x)= {x =-24 > g,(1)= 0 
K =-34 > g,(1)=0 
k=-2>g()=0  g'()=0 par récurrence Vj>0 g;'(1)=0= y(z)=z 


= Vi>3 g;(1)=0 


(D = _À 
K=—2À > Ska E ()=64c, = par récurrence Viz) TER EIER 
g(1)=0  g,'(1)=0 
g.(1)=0 g'(1)=1240,=4 
A7 PE 
K =-34 >< g,(1)=0 CAUSE par récurrence Vj>2 & 0-02 19-442) 
32+K 
1)= 0 ')=62°C, =0 
g;(1) g;"(1) ae? 


Exemple illustratif de l'absence de terme logarithmique sur l'équation différentielle de Heun 


Considérons l'équation de Heun : 


ebe eh fat 


z-l z-a z(z-1\z-a 


Condition de Fuchs y+0+e=a+fp+1 


Point singulier z=0 P=% qo 0= 1(r)=r +r(y- 1)= r=0 r=l-y 


Point singulier z=1 p= du zs DU > I(r) r SCH 1)=r=0 r=1-6 


Point singulier z=a p,=E€ 4) =0 > I(r) r "Arie - 1)= r=0 r=l-E 


Point singulier z = = r° +r(1- Po)+ 4 = (r- ar- B)= 0>r=a r=ß 


Et regardons spécifiquement la condition pour laquelle pour le point singulier z=a, le terme 
logarithmique sera absent. Étant donné qu'autour de la singularité z=a l'écart en valeur absolue 
entre les deux racines indicielle est de |1-e], il suffit que € soit un entier négatif égal à -N pour les 
deux racines soit respectivement : rı=1+N et r,=0. 


Les fonctions f; sont définies comme suit : 


> f(r)=r(r-1+6) 


Lorsque e=-N, r.=1+N et r,=0, la condition d'annulation du terme logarithmique s'écrit sous la 
forme : 


AN) HU DNA — Al lä) 
AN) HNH GN — 0) Zell 
0 GIN f(N-2) … $ 


Fa =n =0)= 0 mn Gol CR 
AO L0) 
|o 0 o Aa) A0] 
AN)=-N 
Avec Al. po + DEE ga 


a'(a-1) 
Lorsque £=0 alors F (0)= f(0)=0< q(a-1)+(-q+aB)a=0< q=a of . 


Lorsque e=1 alors p (0)= E Got = f(1)£(0)- n0) (0)=0. Le numérateur simplifié de cette 
AQ) A0) 
expression donne une équation du second degré en q : 
qg -q{i-y+a(y+8-2+2aB))+a aB(a(y+5-1+aB)-y)=0 
Etant donné que la relation de Fuchs a court : Y + ô-2=a+P Cette équation devient : 


qg -q{i-7+a(a+B+2af))+a apla (1+aX1+B)-7)=0 


AC) Gu) A0) 
Pour ed alors F (r=0)=| f (2) 


© 
> > 
= 
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il s'agit de trouver les racines d'un polynôme du troisième degré en q, sachant que la relation de 
Fuchs a court: Y +8 -3=0+P , L'expression de polynôme en q° est la suivante : 


d'(B(1+B)+2a(1+58+38°)+a°(2+68+3p°)+ 
q-q 83y -4-a (1+3(a +8 +ap)))+a| +a(a(4+48 -4y -6p 7)-4B(y -1))+ 

+2(y Us -2) 
-a apla’ (a? +3a+2\B°+38+2)-ay(4+4a +46 +3aB)+2y(y -1))=0 


Cela coïncide avec l'expression (2.7) obtenue par K.Takemura dans un article de 2022 « Heun’s 
equation, generalized hypergeometric function and exceptional Jacobi polynomial ». Cette 
expression peut aussi s'écrire sous la forme (voir article de 2018 T.A.Ishkhanyan, T.A.Shahverdyan, 
A.M.Ilshkhanyan « Expansions of the Solutions of the General Heun Equation Governed by Two- 
Term Recurrence ») : 


(la-a aB) +(q-a aB\4a-2-a(3+a+pB)+ y)+2a(a-1)ag\q-a of —2a(1+ a + B)+2(y-1))+ 
+2a(a—1(q-a aB\1+aY1+B)=0 


Pour construire les deux développements indépendants pour les deux racines indicielles r,=1+N et 
r,=0 de l'équation de Heun autour de la singularité z=a (ici apparente): 


Date lb Seed) 20-00 EE 
ji WÄER 


Sa) 


Zäit 


sten 
T KN g,"(r) 
0 J=0 


La relation de récurrence pour les coefficients du développement est la suivante : 
jeN g(r) Leiwt e, Jrltlrzuy —1)+...+ a (r)f, fr +1)+ Cf(r)= 0 


Cette expression de la relation de récurrence conduit à un développement de la solution 
suffisamment compliqué pour ne pas poursuivre le calcul plus avant. 


Mais on peut simplifier la récurrence en exprimant l'équation différentielle sous la forme suivante : 


Rite) D Ce af? Cr iser: SEH 


z z—-l z-a 


Soit Po (z)y"(z)+ D (z)y'(z)+ Pa (z)y(2) =0 où Dy (z), D (z) Pa (z) analytique autour de z =a 


e e D H D z de D 
Le développement de la forme : Ga) ce, (z-a conduit à une relation de récurrence à trois 
J 
j=0 


termes, comme suit : 


Pete +R;c =0 
P,=(r+a-1+j\r+B-1+ j) 
O, =-(r+ j} -4-(r+jXr-1+7+e)+al(r+ j} + ap len Aler B-1+6) 


R; =aļa—-1\r+1+jr+e+j) 


E 


L'équation indicielle se retrouve en annulant le degré minimal soit celui du terme (z-a)' : 
Terme (z-a) =0< P, =a(a-1)r(r+e+1)=0 r(r+e-1)=0 


D'autre part pour la racine indicielle r=0, le coefficient R; s'annule lorsque j=-E. Pour que le 
coefficient cı demeure fini et que le développement d'une solution de la forme >. c (z-a) existe 
J 
j=0 
il convient que l'application de la relation de récurrence pour j=-E=N (lorsque £ est un entier l'écart 
des racines indicielles est alors un entier) soit vérifiée : 


E= -N P Cea + Q C; + Ris = 0 


P,=(a-1-EsXB-1+) Q,.=aaB-qg+e(e+1\2a-1)+e(y+e-a(a+B+1+e) R,=0 
S(a-1-sKB-1+ek., Loof -q+e(e+1Ï2a-1)+e(7 +£ -a(a +B+1+6))c., = 0 


Si par construction les deux types de développement SOC eg d(z-a) et 
j=0 


MoS c(z-a) sont indépendants alors la condition précédente garantit qu'il y a absence 
J 

j=0 
complète du terme logarithmique dans la solution indicielle r=0 (racine la plus petite) autour de 
z=a. 


Il faut juste savoir qu'en réalité une solution autour de z=a peut être est construite sous la forme 
d'un développement fini de fonctions hypergéométriques 2F1, ce que nous verrons dans l'exposé 
plus détaillé consacré aux fonctions de Heun. Comme cette solution est probablement 
indépendante de la solution indicielle r=1-e, alors là encore si l'on retombe sur la même condition 
polynomiale conduisant à des valeurs spécifiques de q racine de ce polynôme, on peut affirmer que 
c'est bien la condition d'absence du terme logarithmique. 


